arXiv:1502.02904v3 [cond-mat.supr-con] 13 May 2015 


Journal of the Physical Society of Japan 


DRAFT 
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We study surface states of high-Tc cuprate superconductor YBCO using the bilayer t - J 
model. Calculations based on the Bogoliubov de Gennes method show that a flux phase that 
breaks time-reversal symmetry (7“) may arise near a (110) surface where the ti;t 2 _j, 2 -wave 
superconductivity is strongly suppressed. It is found that the flux phase in which spontaneous 
magnetic fields in two layers have opposite directions may be stabilized in a wide region of 
doping rate, and split peaks in the local density of states appear. Near the surface, spontaneous 
magnetic field may not be observed experimentally, because the contributions from two layers 
essentially cancel out. This may explain the absence of local magnetic field near the (110) 
surface of YBCO, for which the sign of T violation has been detected. 


1. Introduction 

In high-Tc cuprate superconductors, spontaneous violation of time-reversal symmetry 
(T) has been observed in various kinds of experiment. One of the famous example is the 
peak splitting of zero bias conductance in ab-oriented YBCO/insulator/Cu junction. This has 
been interpreted as a consequence of the occurrence of second superconducting (SC) order 
parameter (OP) near the surface, which has symmetry different from that in the bulk.^ For 
this type of surface state, spontaneous current would flow along the surface, and a magnetic 
field should be generated locally. However, experimental evidence for such magnetic fields is 
still controversial.^’®^ 

The present author has studied the (110) surface state of high-Tc cuprates based on the 
Bogoliubov-de Gennes (BdG) method applied to a single-layer t - J model, and found that 
a different kind of T -breaking surface state, flux phase, can occur. The flux phase is a 
mean-field solution to the t-J model in which staggered currents flow and the flux penetrates 
a plaquette in a square lattice.This state has free energy higher than that of the dyi_y 2 - 
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wave SC state except very near half filling, so that it is only a metastable state in uniform 
systems.Near (110) surfaces J;|. 2 _y 2 -wave SC state is strongly suppressed, and then the 
flux phase may arise locally leading to a 7“-breaking surface state. However, the doping region 
in which T violation occurs was much narrower than that observed experimentally in YBCO, 
if we use an effective single-layer model. 

Later reexamination using a bilayer t - J model that describes the electronic states of 
YBCO more accurately have shown that the flux phase may occur as a metastable state in 
a doping region much wider than that for the effective single-layer model. For the bilayer 
t - J model, there may be two types of flux phase in which the directions of the flux in two 
layers are the same or opposite, and a phase transition occurs from the latter to former as the 
doping rate increases. We call the former (latter) one as a type A (B) flux phase. If the type 
B flux phase occurs near the (110) surface, the spontaneous magnetic field should be very 
small, since the contributions from two layers essentially cancel out. This may explain why 
no magnetic field is observed in some experiments for the (110) surface state of YBCO. 

In this paper, we study the (110) surface states of YBCO system that are described by the 
bilayer t - J model. Spatial variations of the OPs are treated using the BdG method,and 
we will show that the flux phase can occur in a wide region of the doping rate when the SC 
order is suppressed. The local density of states (LDOS) is also examined to see whether the 
splitting of the zero-energy peak occurs in agreement with experimental results. 

This paper is organized as follows. In Sect. 2 the model is presented and the BdG equa¬ 
tions are derived. Results of numerical calculations are described in Sect. 3, and Sect. 4 is 
devoted to summary. 


2. Bogoliubov de Gennes Equations 

We consider the bilayer t - J model on a square lattice whose Hamiltonian is given by 
// = //i + //2 + with 


= 


y/.o- (j,t) 


( 1 ) 

( 2 ) 


y'/.o- i 

where the transfer integrals (in plane) tjt are finite for the first- (t), second- it'), and third- 
nearest-neighbor bonds {t"), or zero otherwise. J (7x) is the inplane (interplane) antiferro¬ 
magnetic superexchange interaction, and {],£) denotes nearest-neighbor bonds.The in¬ 
terplane transfer integrals are chosen to reproduce the dispersion in k space,^°^ = 
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(cos kx - cos kyf-, namely, ”on-site” (?q ), second- = -?o/2), and third-nearest-nearest- 
neighbor bonds (t^ = Iq /4) are taken into account. 

c® is the electron operator for the z-th plane in Fock space without double occupancy, 
and we treat this condition using the slave-boson method'^’by writing c® = 
under the local constraint = 1 at every j site. Here (bj) is a fermion 

(boson) operator that carries spin cr (charge e); the fermions (bosons) are frequently referred 
to as spinons (holons). The spin operator is expressed as S® = \ 'Lajs fjT■ 

We decouple the Hamiltonian in the following manner.The bond OPs in plane 
{bf^bf) are introduced, and we denote for nearest- 

neighbor bonds. The interlayer bond OP is defined as x] = ff^)- Although the 

bosons are not condensed in purely two-dimensional systems at finite temperature (T), they 
are almost condensed at a low T and for finite carrier doping (d > 0.05). Since we are inter¬ 
ested in the low temperature region {T < 10~^7 ~ lOK), we treat holons as Bose condensed. 
Hence, we approximate {Af) ~ ^^6 and {A'l’^bf) ~ d (d being the doping rate), and replace 
the local constraint with a global one, 2 Yj,Afj^ fjl) = 1 - d, where N is the total number of 
lattice sites within a plane. This procedure amounts to renormalizing the transfer integrals by 
multiplying d, e.g., t td, etc., and rewriting as In a qualitative sense, this approach 
is equivalent to the renormalized mean-field (MF) theory of Zhang et (Gutzwiller ap¬ 
proximation). The spin-singlet resonating-valence-bond (RVB) OP on the bond {j, 1) is given 
as A® = (iJ;V®-/;f//;V2. The interlayer RVB OP is defined as A| ^ 

Under the assumption of the Bose condensation of holons, Ajj is equivalent to the SCOP. 

We treat a system with a (110) surface, and denote the direction perpendicular (parallel) 
to the (110) surface as .r (y). The .r coordinate is given as xj = jxa where a = a'/ V2 with 
a' being the lattice constant of the square lattice. In order to describe the Flux phase and the 
SC state, x^A^ = and = A*.'^., _ „ are defined for the z-th plane. We assume that 

the system is uniform along the y direction, and consider the spatial variations of OPs only in 
the X direction. By imposing the periodic boundary condition for the y direction, the Fourier 
transformation for the y coordinate is performed.^^^^^^ (Hereafter we write jx simply as j, and 
take a = 1.) Then the MF Hamiltonian is written as follows 

^MF = EE ^](k)hji(k)^i(k), (3) 

k j,t 

with 'Pj(k) = (fj^^\k), fj^^\-k), fj^^\k), fj^\-k)), and k is the wave number along the y di- 
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rection. The matrix hij{k) is given as 

4"(t) f„ 

f; - e ;® 

fj, F'fik) 

fj, , 

where 

= -dpiju + 2t'6 cos 2k) 

-6jj-i[2t6cosk+ 

-6jj,,[2t6cosk+ + ix^Te'^)] 

+ 5jj+2)(t' + 2.t" eos 2k)6, 

6ji(k) = -6ji[(to + 2^2 cos 2k)6 + 

-(Sjj -2 + Sjj+ 2 )(t 2 + 2t^ eos 2k)6, 

F^(I0= + 



fji- 

with n being the ehemieal potential, 

We diagonalize the MF Hamiltonian by solving the following BdG equation for each k, 

^ hji(k)uin(k) = En(k)Ujn(k), ( 6 ) 

i 

where En(k) and Uj„{k) are the energy eigenvalue and the corresponding eigenfunction, re¬ 
spectively, for each k. The unitary transformation using Ujn(k) diagonalizes the Hamiltonian 
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'Hmf, and the OPs and the spinon number at the j site for the layer 1 can be obtained as, 
{nf) = 77 y [|M4;-3,n(^)|V(^n(^)) 

k,n 

+ \u4j-2,n{k)f[^ - f{En{k))\ 

^ Kn (7) 

+U,j^2,n{k)Kj_2,nik)e^'\^ " /(^nC^)))], 

■'’ k,n 

/E tytl \ n 

+Kj_ 2 ^^{k)uAj+i,„{k)e-'’^] tanh (^^)]> 

where Nx (Ny) and / are the number of lattice sites along the x (y) direction within a plane, 
and the Fermi distribution function, respectively. The OPs and the spinon number for the 
layer 2 are obtained by replacing the subscripts, (4j - 3) ^ (4 j - 1), (4 j -2) ^ (4 j), etc., in 
Eq. (7). The interlayer OPs are given as, 

y " E [ulj-2,n(E)UAj-l,n(E)fiEnik)) 

k,n 

+U4J-2Ak)ulj,^(k)(l - f(E„(k)))], 

1 r 1 (^) 

A| = — 2^ [u4j-3Ak)KjJk) + ulj_2/k)U4j-Ln(k)\ 

^ k,n 

xtanh(^). 

V 2T > 

The d- and 5-wave SCOPs are obtained by combining A*^'*^s: - A^”^ -l-A^^j - 

A^>)/4 and A!p(j) = (Af ^ + Af > + A^ + A^;)/4. 

3. Surface States and Local Density of States 

In this section we present the results of numerical calculations for surface states. The 
procedure of numerical calculations is the following. We diagonalize the Hamitotonian Hmf 
with the OPs substituted in matrix elements, and the resulting eigenvalues and eigenfunctions 
are used to recalculate the OPs. This procedure is iterated until the convergence is reached. 
For the system size, = 200 and Ny = 100 are used throughout. The band parameters are 
chosen after Ref. 30; tjJ = 2.5, t'jt = -0.3, t"jt = 0.15, = 0.15, and J^jJ = 0.1. These 

parameters were chosen to reproduce experimental results for YBCO.^°^ We restrict ourselves 
to the case of low temperature, T = 10“^7 (~ IK). 

The spatial variations of the OPs for d = 0.15 are shown in Figs. 1 and 2. It is seen that 
the J-wave SCOP is suppressed near a (110) surface. The imaginary parts of the bond OPs 
(Imy^) are finite there, and Im;^' for different layers have opposite signs. This means that the 
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flux phase arises leading to a 7“-breaking surfaee state. Spontaneous eurrent flowing along 
the surfaee is given as, 




( 9 ) 


with 00 = h/2e being the flux quantum. (In prineiple, there is a term proportional to the veetor 
potential in Jy, but we negleet it for simplicity.) From this equation, we see that the directions 
of the currents and those of the flux in two layers are opposite (type B flux phase). In this case, 
the spontaneous magnetic field near the surface will be very small, since the contributions 
from two layers essentially cancel out. Then it will be hard to observe it experimentally. 
Small imaginary parts of the 5-wave SCOP (Im A^) also appear near the surface, and their 
signs in two layers are also opposite. 



Fig. 1. (Color online) Spatial variations of the SCOPs for 6 - 0.15. Here x is measured in units of lattice 
spacing a. Note that all OPs are nondimensional. 



X 


Fig. 2. (Color online) Spatial variations of the bond OPs for 6 - 0.15. 


The results for other values of 6 show qualitatively the same behavior; the absolute values 
of Im^ are larger (smaller) for smaller (larger) 6, and the surface flux state persists to d ~ 0.3. 
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In MF calculations for uniform systems, the type B flux phase arises for 6 < 0.15, and the 
transition to the type A state occurs as 6 increases. The latter state persists to 6 = 6c ~ 0.2 
in uniform systems.On the contrary, in the present BdG calculation, only the type B phase 
occurs, and 6c is much larger (5^ ~ 0-3) than that in the uniform case. This is because the 
incommensurate flux phase, which is not taken into account in the uniform case, may be 
possible in nonuniform cases, and the type B incommensurate flux state has free energy lower 
than that of the type A incommensurate state. 

For larger Jj^, xj may be a complex number.However, 7^ m that case should be unre¬ 
alistically large, and xj is real for the parameters appropriate for YBCO. 

In BdG calculations, the type A flux phase may be obtained as a metastable state that has 
free energy higher than that of the type B state. In this state, ImA^ in two layers have the same 
sign, in contrast to the case of type B phase. This indicates that ImA^ is induced by Im;^f, and 
its sign is determined by the latter. In the type A state, the imaginary part of the interlayer 
pairing OP, ImA-^, is finite. Since A-^ has the same symmetry as the inplane i’-wave SCOP, 
A^‘\ there is a bilinear coupling term in Ginzburg-Landau free energy, 7 (A^/^ -l- Af^)A-^, with 
7 being a coupling coefficient. This induces ImA-^ once ImA^'^ becomes finite. In the type B 
phase, however, A^ vanishes because A^/^ = -Af\ 

Next we study the LDOS. The LDOS at the j site of the layer 1 is given as 


N:(j, E)= {\uAj-^,n{kt5{E - EM) 


k,n 


( 10 ) 


-t-\UAj-2,n{k)\^6{E + En{k))^, 

and the LDOS for the layer 2 is obtained by replacing the subscripts, (4j - 3) ^ (4j - 1), 
(4 j -2) (4 j). In numerical calculations we replace the 6 function by a Lorentzian with the 

width 0.017. In Figs. 3-5, the LDOS at the surface and in the bulk are shown for 6 = 0.10, 
0.15, and 0.20. (The LDOS for the layer 1 and 2 are the same.) It is found that the splitting 
of peaks occurs in agreement with the experiment. The height of the peaks become larger 
when 6 gets smaller, while the peak splitting, AE, changes only slightly in a nonmonotonic 
way; AE = 0.07637,0.09037, and 0.0111 J for 6 = 0.10,0.15, and 0.20, respectively. 

In order to understand the physical origin of the peak splitting in this model, we show the 
LDOS at the surface of a state with only 7-wave SC order, and that with a surface flux phase 
as well as bulk 7-wave SC order (i.e., without 5 -wave SCOP) in Fig.6. Here the parameters 
are the same as those used in Fig.3. It is seen that the peak splitting occurs as long as the flux 
phase is present. This indicates that the flux-phase order, not the second SCOP ImA^, is the 
necessary ingredient for the peak splitting. We note that it is not possible to have a state with 
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Fig. 3. (Color online) LDOS at the surface and in the bulk for 5 = 0.10. 



Fig. 4. (Color online) LDOS at the surface and in the bulk for 5 = 0.15. 



Fig. 5. (Color online) LDOS at the surface and in the bulk for 5 = 0.20. 


an 5-wave SCOP without the flux phase in the present model. Next we show the LDOS of a 
state with bulk (metastable) flux-phase order (type B) in Fig.7. All SCOPs are set to be zero, 
and the parameters are the same as in Fig.3. The LDOS in the bulk has broad peaks, and one 
of the peaks shifts near to £ = 0 at the surfaee. By eomparing Figs. 6 and 7 with Fig.3, we 
ean see that the peak strueture of the latter is mainly due to the flux phase, and the J-wave SC 
order also contributes to the behavior of the LDOS. 
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Fig. 6. (Color online) LDOS at the surface of a state without s-wave SCOP for 5 = 0.10. 



Fig. 7. (Color online) LDOS at the surface and in the bulk for a flux phase with 6 - 0.10. Here all SCOPs are 
set to be zero 


4. Summary 

We have studied the states near the (110) surfaee of high-Tc euprate YBCO that are de¬ 
scribed by the bilayer t-J model. Near the surface, superconductivity is strongly suppressed, 
and the flux phase in which the directions of the flux in two layers are opposite may oc¬ 
cur in a wide doping region. Then T symmetry is violated and the LDOS at the surface has 
split peaks consistent with experimental findings.'^ The spontaneous magnetic field that could 
arise near the surface with T violation will be very small, because the contributions from two 
layers essentially cancel out each other. These results may explain why no magnetic field is 
observed in some experiments for (110) surfaces of YBCO, for which the sign of T violation 
is detected. 
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